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We present the first experimental observation of resonance-assisted tunneling, a wave phe¬ 
nomenon, where regular-to-chaotic tunneling is strongly enhanced by the presence of a classical 
nonlinear resonance chain. For this we use a microwave cavity made of oxygen free copper with the 
shape of a desymmetrized cosine billiard designed with a large nonlinear resonance chain in the regu¬ 
lar region. It is opened in a region, where only chaotic dynamics takes place, such that the tunneling 
rate of a regular mode to the chaotic region increases the line width of the mode. Resonance-assisted 
tunneling is demonstrated by (i) a parametric variation and (ii) the characteristic plateau and peak 
structure towards the semiclassical limit. 

PACS numbers: 03.65.Sq, 42.55.Sa, 03.65.Xp, 05.45.Mt 


Tunneling describes the possibility of a quantum par¬ 
ticle to transmit through a barrier into a region of space, 
which is inaccessible for a corresponding classical parti¬ 
cle. It is a general wave phenomenon. Dynamical tun¬ 
neling describes the tunneling of waves between classi¬ 
cally disjoint regions of phase space, even without an 
energy barrier being present [1]. It occurs in several 
variants [2], e.g., from a regular region to the chaotic 
region [3-9], from a regular region via the chaotic re¬ 
gion to another regular region [10-16], or between two 
chaotic regions [17-19]. It is essential for applications 
in atomic and molecular physics [20-22] , ultracold atoms 
[15, 16, 23], optical cavities [24-28], quantum wells [29], 
and microwave resonators [7, 14]. 

We consider regular-to-chaotic tunneling, where the 
tunneling rate 7 describes the decay jV'reg(t)P oc 
exp(—yt) of a quantum state, initially located within the 
regular region, to the chaotic region. Towards the semi¬ 
classical limit 7 is determined by two main effects: For 
small wave numbers direct regular-to-chaotic tunneling 
typically leads to an exponential decrease of 7 with in¬ 
creasing wave number [3, 6, 30]. For larger wave numbers 
resonance-assisted tunneling (RAT) drastically enhances 
the tunneling rates, causing the characteristic plateau 
and peak structures [31-34]. RAT occurs due to nonlin¬ 
ear resonance chains inside a regular region, see Fig. 1(b) 
(orange lines), which arise due to the Poincare-Birkhoff 
theorem. A combined prediction for direct regular-to- 
chaotic tunneling and RAT is given in Ref. [34]. 

Several experiments were performed demonstrating 
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FIG. 1. (color online), (a) Photograph of the cavity includ¬ 
ing a movable (white arrow) half disk, absorbers (gray), and 
antenna position (red), (b) Poincare section along the upper 
boundary of the classical dynamics of the closed billiard with 
regular tori (red), a chaotic orbit (blue), and a 3:1 nonlin¬ 
ear resonance chain (orange). The gray region indicates the 
absorbing part, (c) Frequency n and (d), (e) width F of the 
regular modes (n,m) = (1,10) (black pluses) and (4,9) (red 
diamonds) vs. half disk position d. Solid orange lines corre¬ 
spond to the 3x3 matrix model, Eq. (3). Insets visualize 
numerically determined regular modes of the closed billiard. 
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direct regular-to-chaotic tunneling [7, 14-16, 23, 26- 
29]. Recently, the coupling matrix element between two 
modes coupled by a nonlinear resonance chain was very 
nicely observed experimentally in the near integrable 
regime [35] in microcavities [27]. The experimental obser¬ 
vation of the enhancement of regular-to-chaotic tunneling 
rates due to RAT, however, has remained open. 

Here we measure the enhanced line width of regular 
modes due to RAT in an open microwave cavity; see 
Fig. 1(a). Couplings between different regular modes 
are caused by a large 3:1 nonlinear resonance chain; see 
Fig. 1(b). We show the influence of RAT in two ways: 
(i) We induce RAT parametrically by a crossing of the 
frequencies of two regular modes under variation of the 
position of a half disk inset, see Figs. l(c)-l(e); (ii) we 
explore the dependence of the width of regular modes 
towards the semiclassical limit, i.e., for increasing fre¬ 
quency. The characteristic plateau and peak structure of 
RAT is observed showing a good qualitative agreement 
with numerically determined tunneling rates using the 
closed system; see Fig. 2. 

For the experimental realization we use a desym- 
metrized cosine [36, 37] microwave resonator [38-42], 
with a large 3:1 nonlinear resonance chain; see Fig. 1(a). 
The area enclosed by the resonance chain is approxi¬ 
mately half of the area of the regular island, such that 
the characteristic plateau and peak structure of RAT is 
observable on top of the experimental line width. Fig¬ 
ure 1(b) shows a Poincare section with perpendicular 
momentum [pyl/p < 1 vs. arclength along the upper 
boundary of the cosine billiard consisting of regular tori, 
a chaotic orbit, and the 3:1 nonlinear resonance chain. In 
the resonator (height h = 1 cm, length I = 60 cm, width 
left wi = 30.96 cm, and right Wr = 24 cm) we placed 
a half disk (radius r = 9 cm) at the bottom boundary 
with an adjustable distance d from the left boundary, see 
Fig. 1(a). It destroys the hierarchical island structure 
and partial barriers in the chaotic region leading to a 
sharp separation of regular and chaotic dynamics. This 
reduces strong fluctuations of the tunneling rates, which 
would hinder the observation of RAT. Chaotic trajec¬ 
tories spread over the whole billiard also hitting the half 
disk, whereas the regular motion is confined to the region 
left of the half disk. The semiclassical eigenfunction hy¬ 
pothesis [43-45] predicts regular and chaotic eigenmodes 
of the billiard. Numerical examples of regular modes are 
shown in Figs. 1(d) and 1(e) (insets). In phase space 
they localize on quantizing tori at the starting point of 
the arrows in Fig. 1(b), as can be visualized by a Husimi 
representation [37]. 

We open the system by introducing a broad band foam 
absorber on the right-hand side [grey shaded region in 
Fig. 1(b)] thus connecting the chaotic part to the con¬ 
tinuum. The surviving chaotic modes concentrate on the 
remaining part of the chaotic region to the left of the ab¬ 
sorber. This enables the observation of regular-to-chaotic 
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FIG. 2. (color online), (a) Frequency difference Azz = 
— 2^(4,m- 1 ) of modcs (l,m) and (4,m — 1) for m = 
8,.. ., 16 and fixed half disk position d = 30.0 cm. The dot¬ 
ted line indicates the value 0. (b) Experimental widths r(i 
(green pluses) of the inner mode for m = 5,..., 16, the square 
root dependence of Fwaii (orange dotted line), and Fantenna 
(blue diamonds) with the dashed line showing the assumed 
linear dependence at the close encounter, (c) Frat (black 
crosses) showing the initial exponential decay and the char¬ 
acteristic plateau and peak structure of RAT. (d) Numerical 
dimensionless tunneling rates 7{i,m) of the closed system (red 
open circles). 


tunneling by measuring the line width of a regular mode. 
A tedious extraction by analyzing many avoided cross¬ 
ings is not necessary [7]. 

The fixed antenna is placed at the upper left part of 
the billiard [red dot in Fig. 1(a)], guaranteeing a good 
excitation of all regular modes. The complex reflection 
amplitude ^n)^) was measured as a function of the mi¬ 
crowave frequency v. In the isolated resonance regime it 
can be described by the Breit-Wigner formula 
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where the sum runs over all modes with eigenfrequency 
z/fc and eigenwidth F^. Here ak = Re(Afc)'0fc(ri)z/>fe(ri) 
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is the residuum of the resonance, ipki^i) the kth wave 
function at the antenna position fi, and Xk the coupling 
coefficient of the antenna. The real part of Xk describes 
the coupling to the continuum and the imaginary part 
induces a frequency shift Ak ■ We will use the term mode 
for these resonances and reserve further on the term res¬ 
onance to the nonlinear resonance chain. 

The regular modes are labeled by two quantum num¬ 
bers where n counts the number of excitations 

in the horizontal and m in the vertical direction. Their 
width contains different contributions 

= TraT + Twall + r antenna? ( 2 ) 

where Trat is the width enhancement due to regular-to- 
chaotic tunneling in the presence of the absorber. The 
width Twaii is induced by the bottom and top plates 
and the side walls. The width induced by the antenna, 
Tantenna = 2 |a(„_„)|, is given by the residuum 

(i) RAT induced by parameter variation. —We first 
show the signature of RAT from parametric dependen¬ 
cies of two close-by regular modes. The modes (n, m) = 
(1,10) and = (4,9) (Fig. 1 insets) are coupled 

via to the nonlinear 3:1 resonance chain, as n' = n -I- 3 
and m' = m—1. We measure Sii{i') for several half disk 
positions d, fixing all other parameters. We fitted their 
frequency v and width T, giving the parametric depen¬ 
dencies shown in Figs. l(c)-l(e). We observe a crossing 
of the frequencies and an avoided crossing of the widths. 
An increase of the width of the mode (1,10), lying deep 
inside the regular region, is observed once the mode (4, 9) 
is close by, which is a clear signature of RAT. 

A quantitative description of all features displayed in 
Fig. 1 can be obtained in terms of a matrix model. 
We start with a 2 x 2 matrix description of the two 
regular modes with diagonal elements Ei — i^il2 and 
E/i[d) — i^ 4 :l 2 , where Ei and A 4 are the eigenenergies 
of the uncoupled modes, (1,10) and (4,9). As the inner 
mode concentrates away far off the half disk, we assume 
El to be independent of d. Close to the crossing we as¬ 
sume a linear dependence on the frequency axis of E/i{d) 
as is suggested by Fig. 1(c). The widths 71 and 74 are as¬ 
sumed to be independent of d as they are due to antenna 
coupling and wall absorption. An off-diagonal matrix el¬ 
ement Vs:! describes the coupling via the 3:1 nonlinear 
resonance chain. The line width of mode 1, obtained by 
diagonalizing the 2 x 2 matrix, increases at the point of 
degeneracy of the real parts of the two modes, as exhib¬ 
ited in Fig. 1(e). This is a manifestation of the coupling 
of two regular modes via a nonlinear resonance chain. 

For a full analysis of RAT, however, one additionally 
has to consider the tunneling into the chaotic region. If 
the two-mode model were the full truth, the curves in 
Figs. 1(d) and 1(e) should be just mirror images of each 
other because of the invariance of the trace of the Hamil¬ 
ton matrix. Thus the sum of the two line widths should 
be independent of d, which obviously is not the case. 


Thus there must be (at least) one other mode from the 
chaotic region involved. We assume a linear dependence 
of Ech{d) due to the change of the area to the left of 
the absorber causing an average drift in mode energies, 
described by a Taylor expansion up to linear order in d. 
The width 7 ch is assumed to be constant. 

Because of coupling with mode 4 its width dependence 
with d must be essentially the mirror image of Fig. 1(d), 
since the contribution of mode 1 to the width is only 
marginal. Hence the chaotic mode has to become broad 
at the crossing point, while at the same time the width of 
mode 4 becomes small. This is the phenomenology which 
is known in other contexts as resonance trapping [46, 
47]. It is found whenever an effective Hamiltonian of 
the form El = EIq — iWW^ is involved. Here EIq is the 
Hamiltonian of the unperturbed system, and IF is a iV x 
M coupling matrix, where N is the number of modes 
taken into account and M is the number of open channels 
(here TV = 3 and M « 6 is obtained by dividing the 
length of the absorbing edge by A/2). This suggests an 
ansatz of the Hamiltonian in terms of a 3 x 3 matrix. 


H = 


14:1 
0 


14:1 
E4d) - 
— *14ir,4 


0 ^ 

— *14ir,4 


(3) 


This ansatz assumes (a) that the direct tunneling cou¬ 
pling of mode 1 to the chaotic region is negligible for the 
frequency range of Fig. 1 and (/?) that the direct tunnel¬ 
ing coupling of mode 4 to the chaotic region dominantly 
appears in the coupling of mode 4 to the open channels 
encoded in the second row of IF (and much less in the 
coupling of mode 4 to the chaotic mode in EIq). As the 
real parts of modes 1 and 4 cross, with no indication of an 
avoided crossing [Fig. 1(c)], the off-diagonal elements of 
El have to be either real or imaginary, but not complex. 
After diagonalization of the 3x3 Hamiltonian we have to 
compare the complex quantum mechanical eigenenergies 
E = E — i^/2 with the experimentally measured com¬ 
plex electromagnetic eigenvalues v = v — iV/2 hy their 
complex wave number E/ /2M) = {2ttvIc)‘^, yielding 

F; = 1/2 . (27r/c)2; y = 2vY ■ {2tx jcf, (4) 


where we use units F/2M = 1 and F ^ z/. Fitting the ex¬ 
perimental data gives the solid orange lines in Figs. 1(c)- 
1(e) and yields the parameters of the model [48]. In par¬ 
ticular, this gives for the matrix element Vq-i = 2.3 m“ 2 ^ 
which is of the same order as the prediction for the closed 
system, Fsr^cI = 0.51 m “2 [ 49 ]. 

We would like to emphasize that the third chaotic 
mode is an effective description of the coupling of the sys¬ 
tem to the environment, more precisely speaking to the 
absorber. For the example of Fig. 1 we chose a situation 
where the effective broad chaotic mode crosses almost 
at the same d value where the two regular modes cross. 
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For other quantum numbers m the situation was com¬ 
parable but the crossing of the chaotic state was further 
apart. We find a similar value for Vs-a using a reduced 
2x2 model, as this value is defined mainly by the width 
increase of Fji jq) close to d = 29.8 cm. We considered 
either a constant width r4 using only the range 29.5 cm 
< d < 30.0 cm, or an adjusted cosine or Gaussian width 
dependence r4(d) for the whole range of d. This shows 
that the signature of RAT in Fig. 1(e) can be quantita¬ 
tively understood. 

(ii) RAT plateau and peak structure. —A key feature of 
RAT is the plateau and peak structure on top of the ex¬ 
ponential decay of the tunneling rates towards the semi- 
classical limit. We observe this signature of RAT with 
fixed half disk position d = 30.0 cm. The plateau should 
appear when in addition to the inner mode (l,?7i) the 
outer mode (4,m — 1) starts to exist (m > 8). A peak 
is expected at the crossing of these modes for increas¬ 
ing TO, as seen in Fig. 2(a) for the frequency difference 
Ah' = ~ ^(4,m-i) around 6 GHz. There the cor¬ 

responding quantizing tori in phase space are symmet¬ 
ric with respect to the nonlinear resonance chain. From 
the measured widths [green pluses in Fig. 2(b)] 

for TO = 5,. .., 16 we subtract, according to Eq. (2), the 
width contributions F^an and Fantenna (discussed below) 
and the resulting contribution Frat is shown as black 
crosses in Fig. 2(c). For low frequencies up to 3.5GHz 
an exponential dependence is obtained, as expected from 
direct regular-to-chaotic tunneling. RAT becomes visible 
from 4 GHz on as a broad plateau (to > 8) ending in a 
peak structure around 6.5 GHz. 

For numerical comparison, tunneling rates were deter¬ 
mined using the closed system in a similar way as de¬ 
scribed in Ref. [7]. Extending the chaotic region by vary¬ 
ing the height of a rectangle that was attached to the 
bottom of the cosine billiard at x > 30 cm, tunneling 
rates 7(i,m) are determined by evaluating avoided cross¬ 
ings between the regular mode (1 ,to) and chaotic modes 
[Eig. 2(d)]. Qualitatively these rates are in good agree¬ 
ment with the experimental results Frat in Fig. 2(c). 
A quantitative comparison is difficult, as Frat depends 
also on the wall and antenna contributions 74 and 7ch 
which differ for each to. However, the experimental and 
numerical frequency ranges for the initial exponential de¬ 
cay, the plateau, and the peak are identical. 

Let us finally discuss the wall and antenna contribu¬ 
tions Fwaii and Fantenna in Fig. 2(b). The wall absorp¬ 
tion Fwaii has a square root frequency dependence due to 
nonperfect conductance of the cavity walls [50]. It was 
extracted from measurements of the closed cavity, where 
the absorber was removed and the right part of the bil¬ 
liard was closed. We find Fwaii = 0.30 MHz^/j^/GHz; 
see the orange dotted line in Fig. 2(b). The width 
hantenna = induced by the antenna is shown 

as blue diamonds in Fig. 2(b) and is extracted by fit¬ 
ting a Lorentzian to the measured resonance. Around 


6 GHz, where the frequency difference Az/ in Fig. 2(a) 
is small, modes 1 and 4 are coupled, such that the ex¬ 
tracted Fantenna IS affected. In this regime we therefore 
interpolate it by a linear increase; see the dashed line in 
Fig. 2(b). The overall monotonic increase of Fantenna is 
due to two effects: (a) the antenna length is small com¬ 
pared to the wavelength leading to small and increasing 
antenna coupling Re(A), and (/3) the distance of the an¬ 
tenna to the wall is smaller than a quarter of the wave¬ 
length such that the wave function increases monotoni- 
cally at the antenna position. 

In this Letter, we experimentally observe RAT in a 
generic mixed phase space using an opened microwave 
billiard. We demonstrate RAT by the width increase of 
the exemplary mode (1,10) when it is crossed by mode 
(4, 9). With a 3 X 3 matrix model the coupling matrix 
element Van is extracted, which is of the same order as 
the theoretical value. Second, measuring the widths of 
the modes (1 ,to) for increasing to reveals the regimes of 
direct regular-to-chaotic tunneling with an exponential 
decay and RAT with its characteristic plateau and peak 
structure. This experimental work motivates future the¬ 
oretical studies for a better quantitative description of 
RAT, in particular using semiclassical methods [4, 9]. 
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